NASA TECHNICAL 
MEMORANDUM 


NASA TM X- 64539 


DESIGN AND DEVELOPMENT OF A 
LUNAR RADIOMETER 


By Walter Fountain 
Space Sciences Laboratory 


July 28, 1970 


NASA 


CASE FI 

Copy 



George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


MSFC - Form 3190 (September 1968 ) 



TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

TM X -6453 9 . ..... .... 

3. RECIPIENT’S CATALOG NO. j 

4. title and subtitle 

Design and Development of a Lunar Radiometer 

5. REPORT DATE 

July 28 1 970 

6. PERFORMING ORGANIZATION CODE 

7. AUTHOR(S) 

Walter Fountain 

8. PERFORMING ORGANIZATION REPORT # 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

10. WORK UNIT NO. 

1 1. CONTRACT OR GRANT NO. 

13. TYPE OF REPORT & PERIOD COVERED 

Technical Memorandum 

j 

12. SPONSORING AGENCY NAME AND ADDRESS 

14. SPONSORING AGENCY CODE j 

15. SUPPLEMENTARY NOTES 

Prepared by Space Sciences Laboratory, j 

Science and Engineering Directorate . I 

16. ABSTRACT ! 

I 

This report describes the design and presents a performance analysis of an infrared 


radiometer built by the Space Sciences Laboratory. The infrared radiometer is used with an 
astronomical telescope in the study of thermal emission from the moon. Calibration and data 
reduction equations are given, and an error analysis performed for data obtained with the 
system is presented. 


! 


17? KEY WORDS 


18. DISTRIBUTION STATEMENT 

Infrared radiometer 



j 

Astronomical telescope 



! 

Thermal emission 


ANNOUNCE IN STAR. ! 

Calibration and data reduction 



! 

Error analysis 
Lunar temperatures 


i, 3 1 

/y t* *'•" y i 

19. SECURITY CLASSIF. (of this r®pcm 

Unclassified 

20. SECURITY CLASSIF. (of this pag®) 

Unclassified 

21. NO. OF PAGES j 22. PRICE 

! 

44 | $3.00 


MSFC - Form 3292 (May 1969) 


















TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

RADIOMETER DESIGN 2 

Detector Selection and Optical Design 2 

General Performance Analysis 7 

Signal Processing 14 

Mechanical Design 14 

CALIBRATION AND DATA REDUCTION 15 

Calibration 15 

Data Reduction 21 

ERROR ANALYSIS 28 

Propagation of Error Equation . 28 

Error Derivation 31 


REFERENCES 


35 



LIST OF ILLUSTRATIONS 


Figure Title Page 

1. Spectral D* for Ge:Hg, Ge:Cu, and Hg:Cd:Te 

detectors . - ....... 3 

2. Atmospheric transmission, 1 through 4ft jam 5 

3. Schematic of detector/dewar 6 

4. Optical path of radiometer 7 

5. Elements of synchronous system 14 

6. Calibration blackbody slide apparatus 17 

7. Relative response of system 19 

8. Spectral radiance of clear night sky 22 

9. Analog plot of radiometric signal from the crater 

Bruce .' 25 

10. Analog plot of lunar drift scan 26 

11. Expanded scale scan of a portion of lunar disc 29 


LIST OF TABLES 


Table Title Page 

1. System Parameters 10 

2. System N(T) Values 23 

3. Fractional Error Values 


32 



LIST OF SYMBOLS 


Symbol 


A 

e 



D. 


D * 

A 


d 


E' 


F 

n 

N X (X,T) 

N(T) 


N(T ) 
m 

N T (eff) 


N T '(eff) 


NEP, 


Definition 

area of calibration source aperture 

area of detector 

area of telescope objective 

diameter of telescope objective 

spectral detectivity of detector at wavelength A 

distance between calibration blackbody and detector 

image intensity of optical system 

f 

focal ratio of optical system, F n = — 
focal length of telescope 

spectral radiance of blackbody at temperature T 

effective target radiance — integral of the target spectral 
radiance weighted by the system relative spectral response; 
target is assumed to be a black or gray body 

effective radiance of moon at temperature 

effective target radiance — integral of the target spectral 
radiance weighted by the detector relative spectral response 

system noise equivalent radiance — that value of effective 
target radiance which gives rise to a signal -to -noise ratio 
of 1 

noise equivalent power of detector at wavelength A 
absolute power incident on detector 


v 



LIST OF SYMBOLS (Continued) 


Symbol 

P d (eff) 


Rc 


R A (X) 

RR 

A 



s 

S/N 

T 


T 

min 


V 


BG 


V 

Ci 

V 

°2 



V. 


Definition 

power incident on detector weighted by detector's relative 
spectral response 

calibration constant 

spectral responsivity of a given detector 
relative spectral response of a detector 

spectral response of system normalized to unity at a 
wavelength of peak response 

signal output 

ms signal -to -noise ratio 

absolute temperature in degrees Kelvin 

system minimum detectable temperature 

voltage contribution because of background radiation while 
viewing calibration source i 

voltage contribution because of calibration source i 

voltage contribution because of calibration source 2 

voltage because of radiation received from an internal 
reference source 

(V + V_ _,) , sum of voltage contributions from calibration 

Ci B(j 

source i and associated background radiation 
voltage because of sky radiation 


vi 



LIST OF SYMBOLS (Concluded) 


Symbol 


m 


V 


m+s 


a 


AT 


min 


AV 

c 

Af 

e(X) 


T a U) 


T f (A) 


T (A) 
w 


Definition 

voltage output because of the moon alone, (V - V ) 

m+s m 

sum of voltage contributions from the moon and sky 

Fourier coefficient of chop waveform fundamental 

system minimum detectable temperature difference 

rms voltage difference between calibration blackbodies, 
(Vi-V 2 ) 

electrical bandwidth 
spectral emissivity 

average emissivity over a given spectral band 
reflectivity of optical element 
reflectivity of chopper mirror 

product of reflectivity of primary and secondary mirrors of 
Cassegrainian optical system 

spectral atmospheric transmission 

average atmospheric transmission over a given spectral 
band 

actual filter transmission 

average filter transmission over a given spectral band 
actual transmission of dewar window (usually KRS--5) 



TECHNICAL MEMORANDUM X-64539 


DES I GN AND DEVELOPMENT OF A 
LUNAR RADIOMETER 

SUMMARY 


A focal plane radiometer for measuring lunar brightness temperature 
was designed and constructed. The system utilizes a Ge:Hg detector, filtered 
spectrally around a narrow band of wavelengths and centered at 1 1 
micrometers (/mi) . The system was designed to use an f/16 astronomical 
telescope as its collecting optics and is expected to detect temperatures on 
the lunar surface as low as 118°K. An error analysis shows that absolute 
temperatures around the subsolar temperature may be measured to within 
an accuracy of ±7. 9°K with the calibration and data reduction techniques 
described. 

Several scans of the lunar disc near the full moon phase using the 
30-cm reflector of the Space Sciences Laboratory are shown. Additional 
lunar measurements and a description of atmospheric transmission measure- 
ments will be presented in a later report. 


INTRODUCTION 


An infrared (IR) radiometer has been designed to make brightness 
temperature measurements of extended sources in the solar system, and it 
may be useful in studying certain infrared stellar objects such as 
R-Monocerotis or the Magellenic clouds also. The instrument was designed 
to utilize an f/16 Cassegrain telescope as its collecting optics. 

The main object of study in the solar system will be the moon. Using 
the 30-cm telescope of the Space Sciences Laboratory, the system angular 
resolution of approximately 21 arc sec will allow measurements from 
selected areas of about 40 km diameter. Using the Marshall Space Flight 
Center 1. 5-meter telescope at Tucson will give an angular resolution of 
4.5 arc sec, which corresponds to a spot of about 9 km on the moon. 
Measurements of infrared radiation from the planets are possible, but are of 
a limited nature because of the relatively large field of view of the system. 



The system makes infrared measurements in a 2— /im bandpass 
centered near 11 /urn. Since the peak emission of a 264° K blackbody occurs 
at approximately 11 /mi, this might imply that the temperatures of the 
objects of interest are expected to be in this range only. However, any 
object with sufficient emission in this region may be a candidate for study; 
for example, the sun. The selection of the ll-/im region as a design 
specification was determined by the desire to minimize the effects of 
reflected solar radiation, to take advantage of the excellent atmospheric 
window in this spectral region, and to facilitate the conversion of radiance 
measurements to brightness temperatures. 

RADIOMETER DESIGN 
Detector Selection and Optical Design 

The first concern in the design of a radiometer is a determination of 
the approximate strength and wavelength band of the signal to be measured. 
These are necessary criteria for the selection of detector, filters, windows, 
and other optical components. A cursory review of the literature [1] reveals 
that the dynamic range of radiation from an area of the lunar surface may 
be approximated by an extended blackbody source at temperatures of 400° K 
for full moon and about 100°K for new moon. This places the spectral 
emission peaks at wavelengths of 7.2 nm and 28. 9 Atm, respectively. This 
indicates that if a single detector is to be used for both lunar daytime and 
nighttime measurements, it would require adequate spectral responsivity 
covering a large wavelength interval. 

Of the two basic types of detectors available, thermal and photo- 
detective, the photodetective type is considered most appropriate to meet 
all the goals of a spatial scanning type of measurement, where system response 
time is an important parameter. The shorter time constant of the 
photodetective element allows greater latitude in experiment design and 
data acquisition techniques . 

Three commercially available photodetective detectors were considered, 
all theoretically capable of background noise limited performance. These 
were the copper-doped germanium (Ge:Cu), mercury-doped germanium 
(Ge:Hg), and mercury -cadmium -telluride (Hg:Cd:Te) detectors. Typical 
spectral response curves for these detectors are shown in Figure 1. The 
Ge:Cu and Ge:Hg detectors become sensitive to IR radiation when cooled to 
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Figure!. Spectral D* for Ge:Hg, Ge:Cu, and Hg:Cd:Te detectors. 
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temperatures below 20° K s while the Hg:Cd:Te becomes sensitive at 77° K. 

Since these detectors require cryogenic cooling, design requirements must 
include considerations of specific dewar configuration, dimensions, and 
coolant hold times. 

Consideration of the atmospheric effects upon radiation that is 
incident upon the earth from extraterrestrial sources is necessary for guidance 
in bandpass filter selection as well as in signal strength calculations . A 
number of investigators have studied and measured the properties of the 
atmosphere as a function of wavelength for various values of pressure, 
humidity, altitude, and haze conditions. A characteristic absorption curve 
for a vertical path through the atmosphere at sea level is shown in Figure 2. 
This plot shows two general transmission bands that fall within the waveband 
of interest; an excellent band between about 10 and 12 jam and another weaker 
band between 16 and 24 nm. The 16- to 24-gm band is structured with strong 
water vapor absorption bands that lose their strength only at extremely high 
altitudes. A strong 0 3 absorption band occurs at approximately 9. 6 jam. 

Using the 400° K and 100°K values and taking into account the 
atmospheric effects, one can approximate with sufficient accuracy the radiant 
power extremes available, within a given spectral interval, to a detector 
placed in the focal plane of a given telescope. The results of this exercise 
combined with consideration of detector capabilities led to the decision to 
select the Ge:Hg detector and to limit measurements to the 10- to 12-jnm band. 
A schematic of the detector/dewar system is shown in Figure 3. 

The insert in Figure 3 shows the removable detector/aperture unit, 
designed to provide a cold shield for an unvignetted detector field of 
approximately 4 deg. The detector field is thus closely matched to the 
convergence angle of an f/ 16 optical system. 

The optical path of the radiometer is shown schematically in Figure 4. 
Radiation from the telescope optics is reflected off the chopper, passes 
through the KRS-5 window and the interference filter, and is focused at 
the detector surface A. When the two-bladed chopper is in the open position, 
the detector views a reference blackbody, the radiation from the telescope 
is focused at position B where an aperture in a viewing and photographic 
mirror defines the field of view of a 1P21 photomultiplier tube. This channel 
is used to aid in visually correlating the lunar features with the infrared signal. 
Transfer optics reimage the telescope image at point C, where either an 
eyepiece or a camera may be mounted. 
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% ATMOSPHERIC TRANSMISSION 














Figure 4. Optical path of radiometer. 


General Performance Analysis 

Having selected the Ge:Hg detector for the 10- to 12-ptm lunar 
measurements, the characteristics of this detector may be used to evaluate 
certain system characteristics, such as minimum detectable temperature 
and minimum detectable temperature difference. 
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The telescope image irradiance E' from a distant extended graybody 
source, viewed through the atmosphere, is given by 


A A2 

E T = ^4 / N X (X,T) • T fl (X) • P Q HX) • e( A) dA , (1) 

t Aj 


where A* and A 2 are the wavelength limits of interest, and N (A, T) is 

A 

the spectral radiance of the source at temperature T. 

The absolute power to a detector placed in the focal plane of a 

telescope is weighted by the spectral transmission of filters, windows, and 

reflectance of any transfer or chopper mirrors in the optical path. These 

factors are included in the following equation for the power on a detector of 

area A . : 
d 

A °° 

p d = ~P ' A d / 1 V X ’ T > ,T a (A) ' p o {X) ' e(A) ,T f (A) ’ T w (A) * P c (A)dA • 

t ° (2) 


The output signal voltage is determined by the spectral distribution of 

the radiant power incident on the detector weighted by the spectral responsivity, 

R (A), of the particular detector used. Thus, 

X 


u 


A j f RAM 
d 0 A 


N x (A.T) 


T a ( A ) 


P 0< A) 


P C ( A) 


T f ( A) 


t w (A) €(A)dA 


(3) 


In many cases, for a preliminary system analysis, it is necessary to use 
published manufacturer's data for detector performance for filter and win- 
dow transmissions, etc. and to use best estimates for such factors as t 

a 

and optical system parameters. The analysis is simplified by rewriting 
equation (3) as 

A, A2 

S ~ P, ( eff) = —%■ * A • T- • t" *t ' p 2 • € f. RR • N (A, T)dA , (4) 

d f / d f w a o ^ X X 

t Aj 
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where the loss factors are averaged over the spectral bandpass of the system 
as defined by the cuton and cutoff wavelengths of the interference filter, and 
RR^ is a characteristic relative spectral response curve of the type of 

detector to be used. 

The expected signal-to-noise ratio for different lunar temperatures may 
be evaluated from equation ( 4) if the detector spectral noise equivalent power 
is known, since for a narrow band system, 

S v s p d (eff) 

N V N . NEP^ ' “ (5) 


Here, A is taken to be the effective wavelength of the narrow band filter 
used with the system, and a is a factor determined by the waveform of the 
chopped signal. 

Equation ( 5) may be expressed in terms of such parameters as 

detector D* , electrical bandwidth, and optical system F -number by 
A 

combining equations (4) and (5) and making the following substitutions: 



and 


NEP = 
A 


A d */ 2 • (Af) 1 / 2 

D* 

A 


This leads to a general system performance equation, 


S/N 


P d (eff) 

NEP, 


7T _ _ _ _ 9 

a = — • T , • T • T * p 
4 f w a o 

• F ~ 1/2 • tt (eff) 
n I 


•p • € • D* • a A • (Af) 1//z 

c T A d 


(6a) 
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where 


_ A 2 

N T (eff)= / RR a • N A (A,T)dA (6b) 

A i 

and is termed the effective target radiance. 

The system performance will be evaluated for target temperatures of 
200° K and 400° K based on the values listed in Table 1. 


TABLE!. SYSTEM PARAMETERS 


T f 

= 

0.55 

f 


0.70 

w 

T 

a 

= 

0. 90 ( 10 to 12 pm) 

V 

= 

p • p • p = ( 0. 97) 3 = 0. 91 
p s c 

D » 

= 

2. 0 x 10 10 (field of view = 4 deg, Af = 1 Hz, A= 11pm) 

V /! 

= 

"*3 / *“2 

( 1. 96 x 10 cm 2 ) 1 ' 2 = 4. 42 x 10 cm (0. 5 mm diameter) 

Af 

= 

1 Hz 

F 


16 

n 

e 

= 

1 

a 


0.45 


N^(eff) has been evaluated numerically for T = 200° K and 400° K and found 
to be 

A2=!2.2^m 2 

N (eff) = f RR • N (A, T 2 oo)dA= 2.03 x 10 4 watts • cm • sr , 

T 200 V-ln 9 1/rvi A A 
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and 


N, 


T 400 


r 2 

(eff) = / RR 

i 7 


N^(X, T4°°)dA - 5. 


07 x 


«3 

10 watts 


“ 2 ; 


-1 


cm 


sr 


Equation ( 6a) may be written as 


S/N= C • N T (eff) , 


where 



r • p 3 • e • D* • A • a f 1 / 2 
a o X d 


• F ’ 2 = 3.80x 10 5 
n 


cm 2 • sr 
watt 


Thus, 


S/N ( 200°K) = (3.80 x 10 5 ) x (2.03 x 10~ 4 ) = 7.7ix 10 1 , 


and 

S/N ( 400° K) = (3.80X 10 5 ) x (5.07x io" 3 ) = i.927x iO 3 . 


Thus, it is seen that lunar surface temperatures of 200° K and 400° K should 
result in signal-to-noise ratios of approximately 77 and 1927, respectively. 

The minimum absolute temperature measurable with the system is 
defined as that temperature which gives rise to an S/N ratio of 1. If equa- 
tion (6a) is set equal to 1 and solved for N^(eff) , the following expression 

for the system noise equivalent radiance, N^,'(eff) , is obtained: 

V (eff) = c" (7) 
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_g 

Us mg the same parameters as above, this is computed to be 2. 63 x 10 watts * 
cm • sr , which corresponds to a radiometric temperature of approxi- 
mately ii8°K. This then would be the T . of the system based on the 

mm J 

assumed conditions. 


Another parameter of interest is the minimum detectable temperature 
difference, AT^, at some given target average temperature, T & . 

Since N T (eff) may be easily evaluated on a computer in 1°K increments, a 

list of ( eff) values versus blackbody temperature may be tabulated over 

any temperature range of interest. This has been done from 100°K to 400° K. 
The change in effective radiance per degree of target temperature change is 
expressed as 

3N (eff) oo 

— / RR X [N(A,T+1) - N(A,T)]dA . (8) 

Now, the ratio 
N '(eff) 

= AT . (9) 

3N r j,(eff) mm 

3T 


and when evaluated for T = 300° K, gives the minimum detectable temperature 
change, AT . , for targets around 300° K in temperature. Thus, 

“6 -2 -i 

AT . TT'(eff) 2.63x 10 watt • cm • sr • °K „ 

mm .= -zr 1 = = 0. 11 K , 

at 300° K 8N ( eff ) 2.50x10 watts • cm • sr 

T 300 


-6 

AT . 2.63 x 10 

mxn = — 

at 200° K 6.5x10 


0. 40°K , 
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and 


AT . 


mm 


at 400° K 


2. 63 x 10 
4. 22 x 10 


0. 06° K . 


The variation in A T , for different temperatures occurs because 
_ min 

N^(eff) is a function of target temperature and radiation wavelength as seen 
from equation ( 6b) . 

It should be noted that this analysis assumes the sky background 
temperature to be about 295° K because NEP^ was measured in a laboratory 

environment. A turbulent atmosphere would give rise to significantly higher 
background noise levels, which would increase the minimum detectable 
temperature as well as AT^. Conversely, a cold, quiescent atmosphere 

would be expected to give improved performance. 

Examination of equation ( 6a) clearly shows the relationship of the 

output S/N to the various system parameters. It is seen that for a given 

target radiance, S/N is directly proportional to detector D* and area and is 

X 

inversely proportional to electrical bandwidth and system F-number. Improve- 
ments in T . and A T . would require a detector having a better D* or 
min min A 

reduced electrical bandwidth. This latter change would, of course, reduce 
system response time and, therefore, data rate. 

It should be mentioned that for a background noise limited detector, 

D* is inversely proportional to the detector solid angular field of view and 
X 

background temperature. If the detector field of view is cold-shielded to 

match the exit slope angle of the collecting optics, then D* effectively be- 

A 

comes a function of system F-number. S/N would not actually vary as the 
square of the F-number in this case, since a decrease in F^ requires an 

increase in the exit slope angle of the optical system. This in turn requires 
a corresponding increase in the detector cold-shielded field of view, which 
results in an increase in background noise. This fact would have to be con- 
sidered in any estimate of S/N improvement to be gained by changing to a faster 
optical system. In many instances, instruments are designed with respect to 
an existing telescope that is available to the experimenter. 
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Signal Processing 


The ac signal from the detector is brought into a phase -lock amplifier 
and synchronously demodulated to give a 0- to ± 10-Vdc output. The output 
is recorded on a strip-chart recorder or magnetic tape. The amplifier has 
a variable post detection time constant, which may be increased for weak 
signals. An electronic zero offset is provided so that small signal differences 
riding on a large basic signal may be amplified and observed. A brief but 
clear explanation of synchronous detection may be found in Reference 2. 

Figure 5 shows a block diagram of the elements of a synchronous detection 
system . 



Figure 5. Elements of synchronous system. 

Mechanical Design 

The following factors have been considered in the mechanical design 
of the radiometer: 

1. It is rigid enough so that bending or flexing will not affect optical 
alignment. 

2. It is sturdy enough to be transported from laboratory to observ- 
atory without affecting optical alignment. 

3. It is lightweight enough to be mounted on a telescope by only two 

people. 
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4. Heat producing elements such as chopper motors are situated 
such that they do not enter the detector field of view. 

5. Light producing elements such as mirror illuminators or reference 
sources are situated so that they do not enter the photomultiplier field of 
view. 


6. The mounting of the dewar is such that the liquid coolant will not 
spill from the open top with the telescope in any position from vertical to 
near horizontal. 

7. The radiometer should be rotatable about the telescope optical 
axis, and tolerances should be such that this rotation will cause negligible 
optical misalignment. This is very difficult to achieve, especially with 
instruments whose mass is not symmetrical about the optical axis . Future 
designs should provide an x-y adjustment of the entire radiometer while 
mounted to the telescope backplate. Under these conditions, any misalign- 
ment caused by instrument sagging or insignificant mechanical tolerance 
may be corrected by vernier screw adjustments. 

8. All element mounts are accessible for adjustment. 

CALIBRATION AND DATA REDUCTION 
Calibration 


The final form of a radiometric calibration technique is usually 
simple, after one has accurately measured the transmission and/or reflectance 
of all optical elements such as windows, mirrors, filters, beam splitters, 
etc. ; has measured the relative response of the particular detector; has 
accurately determined the true field of view (as opposed to a calculated field) ; 
and has determined if changes in detector orientation or ambient pressure 
affect system response. Other effects, such as polarization and saturation 
effects, etc. , should be measured if possible. If extreme caution and dili- 
gence for precision are practiced, one may expect a calibration error on the 
order of i or 2 percent. 

For a radiometer such as this, it is preferable to make measurements 
simultaneously on an entire series of elements as they will be used in the system 
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rather than making individual measurements on separate elements. This is 
because interactions such as interference and multiple reflections can produce 
a significantly different throughput from that inferred by mathematically 
combining individually measured component characteristics. 

The final step in radiometric calibration is a determination of the 
system responsivity, R, defined as the voltage (or current) output per unit 
radiant input. Of several methods available, the one described below was 
chosen for its simplicity and its relative ease of performance under the 
conditions imposed by observatory operation. 

Two temperature -controllable blackbodies of identical design are 
mounted on a mechanical device (Fig. 6) which is placed on the radiometer 
when a calibration is desired. An accurately machined aperture, i. 52 cm in 
diameter, is located in the stationary base of the device and is positioned by 
guide pins directly above the entrance aperture of the radiometer coincident 
with the system optical axis. This aperture serves as the effective calibration 
blackbody aperture and is at an accurately measured distance from the 
detector. The two blackbodies can be moved quickly and accurately over the 
aperture using the slide arrangement shown in Figure 6. After the 
calibration temperatures of the two blackbodies have been set and allowed to 
stabilize for a minimum of 45 minutes, the detector voltage output is 
recorded for each blackbody in rapid succession. This completes the 
measurement portion of the calibration. 

Two blackbodies are used at two preset temperatures, rather than 
setting a single one at several temperatures consecutively. This is because 
significant background radiation is received from the outer regions of the 
detector's field of view, and, as may be seen from equations ( 10) and ( 11) , 
any change in this background radiation during calibration would introduce 
an indeterminable error into the output signal. The relatively longtime 
required for a single blackbody to stabilize at different temperatures is an 
invitation to this type of error. 

The dc voltage output of the phase lock amplifier is derived from the 
radiative contribution of several distinct sources. The peak-to-peak value 
of the square wave ac signal from the detector is proportional to the 
radiation received on alternate half cycles of the optical chopper. On the 
closed portion of the chop cycle, the detector "sees" radiation from the 
calibration blackbody as well as an unknown amount of radiation from the 
surrounding background (in some cases, this background contribution to the 
signal may even be greater than that from the calibration source) . This may 
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Figure 6. Calibration blackbody slide apparatus, 
be symbolized as 


V + 
c i 


V = V 
BG T 


On the open half of the chop cycle, the detector sees the radiation from an 
internal reference source (or in some cases, simply the radiometer cavity 
radiation) which is symbolized by V . The magnitude and phase relation of 

the resulting ac signal is dependent on the difference ( V - V^) , which gives 

zero output when both are equal. 

The rms output from the electronics when viewing calibration source 
(1) at temperature T* is then 


l(V + 
ci 


V BG ) 


- Vj] a 


( 10 ) 
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and when viewing calibration source (2) at temperature T 2a the rms output 
is 


V 2 = 



+ v 


BG 


- Vj]a 


(ii) 


where a is the Fourier coefficient of the chop waveform fundamental and is 
maximum for a square wave. Since the slide arrangement allows the two 
measurements to be made quickly, the assumption is made that the 
contributions to the signal from the background and internal reference remain 
constant during the measurements. Therefore, the voltage difference 
( V| - V 2 ) is a result of the calibration sources temperature difference only. 

Thus, 


AV c = , V, - V,) = 1 ( v ci + V BG - Vj) - ( v 02 + v BG - Vj)] • o. 

^ = l T o,-V •“ ' < 12 > 


where AV is termed the calibration signal. Referring to equation ( 4) , the 
c 

voltage difference between the two blackbodies may be expressed in terms of 
system parameters and radiometric quantities. This leads to 


A • A oo R 

AV c = AP d (eff) • = - d r - R * / -5-^ * ej\) • r„(X) • r,(X) 

P 


'X J n R, 
p 0 X 


t (X) • PJV • [ N ( X, T i) “ N (X, T 2 )]dX . (13) 

W C A A 


The functions r^(X), t w (X), and P Q (X) are all produced by permanent 

elements of the radiometer, and their absolute values are not required as 
long as absolute power measurements are not required. The relative spectral 
response of the system is required, however, to evaluate the integral. This 
function, normalized to its peak value, was obtained with the aid of a Perkin 
Elmer Model E-14 spectrometer and is shown in Figure 7. 
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RELATIVE RESPONSE 



Figure 7. Relative response of system. 
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If t ( A) is taken to be unity over the calibration distance d and 

ci 

e ( A) is assumed to be constant over the spectral interval, these quantities 
c 

may be brought outside the integral. It is convenient to define a quantity 
N(T) , given by 

°° 

N(T) = f - N (A, T)dA , (14) 

0 \ x 

R 

where — — is the system response normalized to unity at the wavelength of 
P 

peak response as shown in Figure 7. Equation (14) has been evaluated 
numerically and tabulated in 1°K increments from 50 to 500° K. 

Rewriting equation ( 13) gives 


A A d 

AV c = -vf — [N(Ti) -N(T 2 )] , (15) 

P 

which forms the basic calibration equation from which the system responsivity 

constant, R , is evaluated. It would be possible to simulate telescope 
A 

P 

reflectance values during calibration by placing two mirrors, coated 
simultaneously with the telescope mirrors, in the calibration optical path, d. 
These mirrors would not be normal to the optical line of sight, however, and 
would require checks against possible polarization effects. 

When viewing extended sources such as the moon, the factor 

occurs in both the calibration and data reduction equations so that equation ( 15) 

A "e 
c c 

may be solved for the quantity A R, . If one lets A R = R and — -rs— = c > 
J 17 d A dAc d 

P P 

then the final form of the calibration equation is given by 


AV 

R = £ 

c c • [N(Ti) - N(T 2 )] 


(16) 
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R has units of volts • cm 2 • watt 1 and is determined by the method described 
c 

in this section for blackbody temperatures and T 2 . A slightly different 

calibration technique would be to set the calibration blackbody to a number of 

different temperatures and plot the corresponding voltage outputs against 

c • N(T) . The slope of the best fit straight line plot would then yield R . 

c 

This method, however, is more time consuming and would be subject to the 
background change errors discussed previously. 

The largest error associated with this method of calibration comes 
from inaccuracies in the measurement of the absolute temperatures of the 
calibration blackbodies . 


Data Reduction 


Two separate measurements are required to extract the moon signal 
from the moon-plus-sky emission signal. As may be seen from Figure 8, 
the spectral radiance of a clear, cool night sky [2] at ll/um is low, well 
below equivalent 0° C blackbody emission even at elevation angles as low as 
10 deg. However, in regions of high humidity and warm temperatures, the 
sky radiance can exceed that of a 0°C blackbody at low elevation angles, as 
seen from the Cocoa Beach plot of Figure 8. The effect of increasing air mass 
on sky radiance and, hence, sky signal voltage may be seen from the curves 
for different elevation angles. Even the presence of unobservable (to the 
eye) clouds or haze can cause a significant variation in sky level readings 
from one part of the sky to another. Any change in general atmospheric 
conditions such as would result from the invasion of a cold or warm front 
or the influx of a humid air mass will show itself as a change in shape on 
plots of sky radiance versus elevation angle as a function of time. If the 
change is rapid, an obvious change in radiometric output could be observed 
at a given elevation angle. These effects can cause quite a problem in the 
determination of nightly extinction, where data have been taken over long 
observation periods (4 to 8 hours) . 

The contribution to the total signal from sky emission may be 
eliminated simply by observing the sky at a given zenith angle in the moon's 
path and recording the system voltage output. As the moon drifts into the 
field of view, its contribution to the signal is added to that of the sky; thus, 
the difference between the two voltage levels is because of the moon alone. 

As long as this process does not take too long (on the order of 2 or 3 minutes) , 
variations in sky emission may be considered negligible, except as a 
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contribution to system noise . Another effect that may be significant during a 
drift scan is the contribution to the signal from radiation scattered into the 
detector as the bright limb drifts near the radiometer field of view. The 
effect of this scattered radiation would be to increase the sky level reading to 
a level slightly higher than the true sky level. This scattering effect is 
greater at large zenith angles, since the density of scattering particles along 
the slant path is greatest. The magnitude of this effect can be measured by 
pointing the radiometer field of view sufficiently ahead of the travel path of 
the moon to assure that scattering effects are not present and monitoring 
an expanded scale output of sky level voltage as the moon approaches. When 
the surface of the moon is fully in the field of view, errors because of 
scattering effects would be minimal, since, statistically, radiation scattered 
in and out of the field would be about equal for an extended source. 

Data reduction consists of determining a blackbody radiance value 
N(T ) for the moon and comparing this value with the tabulated list of N(T) 

values as determined from equation ( 14) . .A list of system N(T) values 
with their corresponding temperatures is shown in Table 2 for temperatures 
near 400° K. 


TABLE 2. SYSTEM N(T) VALUES 


N(T) 

T 

_ 2 -1 

Watts • cm • sr 

°K 

4.658 x 10" 3 

390 

4. 699 x io’ 3 

391 

4. 740 x 10" 3 

392 

4.78ixl0~ 3 

393 

-3 

4. 822 x 10 

394 

4. 864 x 10" 3 

395 

5. 073 x 10” 3 

400 
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The voltage output of the system for the moon is given by the 
expression, 


V 


(m) 



e N(T ) 
m m 


R 

c 


J 


where 


v = V - v . 
( m) s+m s 


Solving this equation for N(T ) and substituting from equation (16) 

for R gives the desired data reduction equation , 
c 


N(T ) = 
m 


V (V m -c-WTlI -NTOl 

AF 2 • r(a) • AV • £ 
t t cm 


(18) 


AV is measured directly for calibration blackbody temperatures Tj and 
c 

T 2 , and the corresponding values for N(Ti) and N(T 2 ) are taken from the 
listing of N(T) values. V m+g and are the system output voltages 

developed when viewing the moon and sky respectively through a telescope of 
effective collecting area A , focal length f, and mirror reflectance p . 

I l 

The atmospheric transmission is determined separately, as a function of 

air mass, by using the "weak-line approximation" [3] model, and a reference 

point on the moon as the extinction source. N(T ) is the calculated effective 

m 

radiance value of that area of the moon filling the system field of view and is 

convertible to a radiometric temperature, given a value for T . If ~e is 

mm 

taken to be 1 , the temperature is then termed brightness temperature . If the 

goniometric emission characteristics of the area observed is Lambertian, 

then the true temperature would be greater by a factor of i/7 , where 7 

t t 

is the true hemispherical emissivity of the surface averaged over the spectral 
bandpass of the system. 

Two methods of data reduction will be illustrated using the data 
reproduced in Figures 9 and 10. Figure 9 represents an analog voltage plot of 
radiation measured from a selected site on the lunar surface, in this case 
the crater Bruce near the disc center. The radiometer/telescope was visually 
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Figure 10. Analog plot of lunar drift scan. 
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pointed at the site by means of the reflex optics and allowed to track the point 

for approximately 5 seconds. The RA drive is disengaged at point A, and 

the moon drifts out of the field of view. This takes about 1 minute, during 

which time it is assumed that atmospheric emission and transmission remain 

relatively constant. As the moon’s limb leaves the field of view, the signal 

drops rapidly to sky level B, which is recorded for an additional 30 seconds. 

The signal from the moon (V - V ) is determined by direct measurement 
& m+s s 

from the strip chart. Time is recorded so that the zenith angle for the site 

Bruce may be determined from the lunar ephemeris and applied to the 

calculation of air mass for which t is determined. 

a 

(V - V ) and t are entered into equation (18) along with the 
m+s s & 

system constants and calibration data, and N(T^) is calculated for a 40 

km diameter area around the crater Bruce. A computer program has been 
written for the Univac 1108 which will print out the temperature corresponding 
to the nearest NfT^) value that it has in core storage. This is a brightness 

temperature that can be corrected for e, should such value ever become 

l 

available for the area under measurement. The system propagation of error 
analysis in the next section shows that with the present system, the probable 
error in a single brightness temperature measurement of this type is 
approximately ±8°K. 

Figure 10 represents the analog output of the system for a drift scan 
of the disc. The same procedure of scanning is used except the telescope 
is driven to some point in the sky ahead of the moon. The right ascension 
drive is disengaged allowing the moon to drift across the' stationary field 
of view. A full disc scan requires about 2 minutes plus an additional 30 
seconds of sky on each end of the scan. One advantage to this method of 
scanning is that air mass remains constant, simplifying t calculations. 

SL 

During the course of the scan, photographs are taken through the radiometer 
optical system to pinpoint the scan on a separate full phase photograph. 
Assuming a 21 arc sec field of view', about 100 contiguous scans would be 
required to completely scan the moon using this procedure. If an average of 
3 minutes is allowed per scan, then a complete scan program of the full disc 
would require 5 hours . 

A temperature scale may be derived for this type of data in the 
following manner. The data reduction equation is solved by inputting a 
range of values for N(T^) and calculating a corresponding list of 
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(V - V ) values. These values are printed out in either i° or 10° 
m+s s 

increments of temperature versus voltage for each scan. The temperature 
scale may then be applied directly to the strip chart voltage trace as shown 
in Figure 10. Figure 11 shows an expanded scale scan of the central portion 
of the disc using the electronic offset feature of the synchronous detection 
system. 


ERROR ANALYSIS 
Propagation of Error Equation 

A radiometric measurement of temperature actually requires the 
individual measurement of a number of other quantities from which the 
temperature value is calculated. To determine the probable error 
associated with this temperature value, one must know how the inaccuracies 
associated with the determinations of these other quantities contribute to the 
inaccuracy of the result. 

The system propagation of error equation derived below shows the 
contribution of the individual component errors to the overall temperature 
error. Where possible, the probable errors associated with the individual 
components were derived by standard statistical means. On items where 
measurements were not feasible, best estimates of error are used. 

Following the method used by Ingrao [4] , the data reduction equation, 
which is an expression for the radiance of the moon as measured by the 
system, is expressed in general form similar to the Stefan-Boltzmann law as 


N(T ) = A T B (19) 

m m 

where A and B may be assumed constants over the range of temperatures 
expected on the sunlit portion of the lunar disc. Differentiation of this 
equation with respect to gives 

dN(T ) = ABT. . B-i d T (20) 

m (m) m 
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Figure 11. Expanded scale scan of a portion of lunar disc 






Replacing A with 


leads to 


. d N(T ) 
1 # m_ 

B ’ N(T ) 
m 


N(T ) 
m 


m 

dN(T 


dN(T 


N(T ) 
m 


If the errors in N(T ) and are small, the following approximation is 
valid [ 5] : 

AT r N(T ) dT 1 A N(T ) 


N(T ) dT 
m . m 

T dN(T ) 

m m 


A N(T ) 
m_ 

N(T ) 
m 


The weighting factor 


N(T ) dT 
m _ m 

T ’ dN(T 


may be evaluated using the 


pretabulated tables of T versus N(T), and the expression for 


AN(T ) 
m 

N(T ) 
m 


may be determined by differentiating equation (18) . This leads to the 
expression, 


’ dN(T ) 
m m 

(A( V -V ) 
m+s s 

V -V 
m+s s 


AA 

c 


Ae 

c 

t o 

if t 


4-i 

< 

<1 


Ae 

m 

A 

c 

+ 

e 

c 

+ Ci 

f t 

+ 

A 

t 

-V 

e 

m 


At 

a 

* o 

Ap t 

• 9 

Ad 


T 

a 

T ^ 

P t 


d_. 

+ 


Vi-V 2 


/ dN\ 

V dT A 


ATj 

N(T 1 )-N(T 2 ) 


AT 2 

N(Tj) -N(T 2 ) 
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Evaluation of this expression with maximum expected fractional error values 
entered for the components of the second bracket would yield the maximum 
fractional error expected for a single temperature measurement. 


By defining the probable error of a single temperature measurement 
of the moon in a manner analogous to the standard deviation (i.e. , the square 
root of the sum of the squares of the individual probable errors) equation (23) 
may be expressed as 



prob 


N(T ) dT 
m . m 

T *dN(T ) 
m m 



+ 


+ 


A(Vl -V j j A 2 + 

Vi - v 2 y 


dN \ ATj 

dT,/^ N(Ti)-N(T 2 ) 



AT; 

N(Ti)"N(T 2 ) 


2 ^ 


/ 


i/2 


(24) 


In this case, the values entered into the second bracket would be the fractional 
probable errors derived for the individual components. 

Error Derivation 

The probable error has been evaluated based on Table 3. The probable 
error is evaluated at two lunar temperatures, T^ = 400°K and 200°K. 

Blackbody calibration temperatures are taken to be Tj = 400°K and T 2 = 300°K. 
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Component 


Fractional Error 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

( 10 ) 
(ID 


A( V -V ) 

m+s s 

V -V 
m+s s 

AA 
c 

A 


± 0.01 


± 0.001 


Ae 


± 0.01 


Af 

t 

£ t 

ttt 

\ 

Ae 

m 

e 

m 

At 

a 


± 0.01 


± 0.01 


±0.03 


±0.03 




± 0.02 


Ad 

d 


±0.0025 


A( Vi~V;>) 
Vi-V 2 


± 0.01 


ATj 


± 0. 25°K 


(12) AT 2 


± 0. 25°K 






For 


T = 400° K, 
m 


N 


(400) 


(400) 
dT 


5. 073 x 10 


-3 


4x 10^ 


i. 268 x 10 


-5 


(J00) _ 2 37x 10 4 
d N 


d N 


dT 


dN 


= 4.219 x 10" 5 


dT 


(400) 


■^22 = 2 . 447 x 10“ 5 


(300) 


N„ 


(400) 


dT 


(400) 


dN 


[= 0.3048 


N_ = 5. 073 x 10 

T ( 400) 

N = 1.706x 10 

^ (300) 


dN 


dT 


(400) 

dN 


dT 


(300) 


1 

[N(400) -N(300)] 

1 

[N(400) -N(300)] 


4. 219 x 10 
3. 367 x 10 


= 1.258 x 10 


2 . 447 x 10 
3. 367 x 10 


7. 98 x io“ 3 




= ( 0. 305) [ ( 0. 01) 2 + (0.001) 2 + ( 0. 01) 2 + 4(0. 01) 2 
+ (0.01) 2 + ( 0. 03) 2 + ( 0. 03) 2 + 4(0. 02) 2 
+ 4(0. 0025) 2 + ( 0.01) 2 
+ (1.258 x iO _2 x 0.25) 2 
+ (7. 98 x iO~ 3 x 0.25) 2 ] 1 / 2 
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AT 


= (0.305) (0.0648) = 0.0197 = ± 1.97% 


(400) /prob 


or ± 7. 9°K . 


The same analysis assuming T^= 200° K results in a probable error of 
±1.01 percent or ± 2°K. This is because of the fact that for the colder 

N(T ) dT 

temperature, the factor — • ^ - — — is reduced from 0. 305 to 

m m' 

0. 156. This estimate is probably somewhat optimistic, however, since for 
the low temperature radiation, greater uncertainties in the measurement of 
atmospheric transmittance and decreased signal-to-noise ratio of the moon 
signal would cause the probable errors of these components to be increased. 
In both cases, the analysis assumes good observing conditions. 


From the error analysis, it is evident that a major source of error 
could be eliminated by calibrating the system through the telescope optics. 


Af, 


AA. 


AP, 


This would eliminate the factors 4 


and 4 


from 


equation (6) . However, with the present Cassegranian systems, a special 
calibration technique such as the Jones [6] method must be employed. This 
method is difficult and time consuming to set up and perform, and is not 
conducive to observatory conditions where night-to-night variations in system 
responsivity require frequent calibrations. For these reasons such an 
approach has not been adopted. 
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